Mon. Not. Ft. Astron. Soc. 000, ITH2T1 (20091 



Printed 28 January 2009 



(MN M£X style file v2.2) 



The Magnetosphere of Oscillating Neutron Stars in 
General Relativity 

Ernazar B. Abdikamalov 1 ' 2 *, Bobomurat J. Ahmedov 2,3 and John C. Miller 1 ' 4 

1 SISSA, International School for Advanced Studies, and INFN-Trieste, Via Beirut 2-4, 34014 Trieste, Italy 

2 Institute of Nuclear Physics, Ulughbek, Tashkent 100214, Uzbekistan 

3 Ulugh Beg Astronomical Institute, Astronomicheskaya 33, Tashkent 100052, Uzbekistan 

^Department of Physics (Astrophysics), University of Oxford, Keble Road, Oxford OX1 3RH, UK 



Accepted <date>. Received <datc>; in original form <date> 



ABSTRACT 

Just as a rotating magnetised neutron star has material pulled away from its surface 
to populate a magnetosphere, a similar process can occur as a result of neutron- 
star pulsations rather than rotation. This is of interest in connection with the overall 
study of neutron star oscillation modes but with a particular focus on the situation for 
magnetars. Following a pre vious Newtonian analys is of the production of a force-free 
magnetosphere in this way Timokhin et al. (2000), we present here a corresponding 



gcncral-rclativistic analysis. We give a derivation of the general relativistic Maxwell 
equations for small-amplitude arbitrary oscillations of a non-rotating neutron star 
with a generic magnetic field and show that these can be solved analytically under the 
assumption of low current density in the magnetosphere. We apply our formalism to 
toroidal oscillations of a neutron star with a dipole magnetic field and find that the 
low current density approximation is valid for at least half of the oscillation modes, 
similarly to the Newtonian case. Using an improved formula for the determination of 
the last closed field line, we calculate the energy losses resulting from toroidal stellar 
oscillations for all of the modes for which the size of the polar cap is small. We find that 
general relativistic effects lead to shrinking of the size of the polar cap and an increase 
in the energy density of the outflowing plasma. These effects act in opposite directions 
but the net result is that the energy loss from the neutron star is significantly smaller 
than suggested by the Newtonian treatment. 

Key words: stars: magnetic field - stars: neutron - stars: oscillations - pulsars: 
general 



INTRODUCTION 



Study of the internal structure of neutron stars (NSs) is of fundamental importance for subatomic physics since these 
objects provide a laboratory for studying the properties of high-density matter under very extreme conditions. In particular, 
there is the intriguing possibility of using NS oscillation modes as a probe for constraining models of the equation of state 
of matter at supranuclear densities. It was suggested long ago th a t if a NS is oscillating, then traces of this might be 
revealed in the radiation which it emits ( Pacini fc Rudermanl 19741 : Tsveanl 19751 : Boriakofl 19761 : Bisnovatvi-Koean 1995 : 
Ding fc Chen3ll997l : [Duncan! 19981 ). Recently, a lot of interest has been focussed on oscillations of magnetized NSs because 
of the discovery of gamma-ray flare activity in Soft Gamma- Ray Repeaters (SGRs) which are thought to be the very highly 
magn etised NSs known as magnetars (for recent review on the SGRs see IWoods fc Thompson! I200I IWatts &i Strohmaverl 
20071 ). The giant flares in these objects are thought to be powered by global reconfigurations of the magnetic field and it 



has been suggested that the giant flares might trigger starquakes and excite global seismic pulsations of the magnetar crust 
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( Thompson fc Duncan! Il995l . bOQll ; ISchwartz et~al1 120051 ; |Puncanlll998l ). Indeed, analyses of the observations of giant flares 
have revealed that the decaying part of the spectrum exhib its a number of quasi-periodic oscillations (QPOs) with frequencies 
in th e range from a few tens of Hz up to a few hundred Hz ( Israel et al. 20051 ; Strohmaver fc Watts!l20o"i ; Watts fc Strohmaver 
20061 ) and there has been a considerable amount of theoretica l effort attempting to identify these with crustal oscillation 
modes i Glampedakis et al.ll2006l : Samuelsson fc Andersson 2007; L cvinll2o'o7 : Sotani et al.ll2007al lbh. While there is substantial 
evidence that the observed SGR QPOs are caused by neutron star pulsations, there i s a great deal of uncertainty about how 
stellar surface motion gets t ranslated into the observed features of the X-ray radiation (|Strohmave rll2008l ; IStrohmaver fc Watts! 
20061 ; iTimokhin et al.l 120071 1. To make progress with this, it is necessary to develop a better understanding of the processes 
occurring in the magnetospheres of oscillating neutron stars. 

Standard pulsars typically have magnetic fields of around 10 12 G while magnetars may have fields of up to 10 14 — 10 15 
G near to the surface. Rotation of a magnetized star generates an electric field: 



SIR 



B 



(1) 



where B is the magnetic field strength, c is the speed of light and SI is the angular velocity of the star with radius R. 
Depending on the rotation velocity and the magnetic field strength, the electric field may be as strong as 10 V cm" 1 and it 
has a longitudinal component (parallel to B) which can be able to pull charged particles away from the s tellar surface, if the 
work function is sufficiently small, and accelerate them up to ultra-relativistic velocities. This result led iGoldreich fc Julian 
( 19691 ) ;o suggest that a rotating NS with a sufficiently strong magnetic field should be surrounded by a magnetosphere 
filled with charge-separated plasma which screens the accelerating electric field and thus hinders further outflow of charged 
particles from the stellar surface. Even if the binding energy of the charged particles is sufficiently high to prevent them being 
pulled out by the electric field, the NS should nevertheles s be surrounded by charged particles produced by plasma generation 
processes (Sturrock 1971 ; Ruderman fc Sutherland 1975h , which again screen the longitudinal component of the electric field. 
These considerations led to the development of a model for pulsar mag netospheres which is frequently called the "st andard 
model" (an in depth discussion and review of this can be found in, e.g.. lMichellll99ll ; iBeskin et al. l ll993l ; lBeskinll2005h . 

Timokhin, Bisnovatyi-Kogan & Spruit (2000) (referred to as TBS from here on) showed that an oscillating magnetized 
NS should also have a magnetosphere filled with charge-separated plasma, even if it is not rotating, since the vacuum electric 
field induced by the oscillations would have a large radial component which can be of the same order as rotationally-induced 
electric fields. One can show this quantitatively by means of the following simple arguments. To order of magnitude, the radial 
component of the vacuum electric field generated by the stellar oscillations is given by 



E° 



—B 



(2) 



where ui is the oscillation frequency and f is the displacement amplitude. Using this together with Eq.JT]), it follows immediately 
that the electric field produced by oscillations will be stronger than the rotationally induced one for sufficiently slowly-rotating 
neutron stars, having 



«*3 



(3) 



For stellar oscillations with £/R ~ 0.001 and u) ~ 1 kHz, the threshold is SI ~ 1 Hz. Within this context, TBS developed 
a formalism extending the basic aspects of the standard pulsar model to the situation for a non-rotating magnetized NS 
undergoing arbitrary oscillations. This formalism was based on the assumption of low current densities in the magnetosphere, 
signifying that the influence of currents outside the NS on electromagnetic processes occurring in the magnetosphere is 
negligibly small compared to that of currents in the stellar interior. This assumption leads to a great simplification of the 
Maxwell equations, which then can be solved analytically. As an application of the formalism, TBS considered toroidal 
oscillations of a NS with a dipole magnetic field, and obtained analytic expressions for the electromagnetic field and charge 
density in the magnetosphere. (Toroidal oscillations are thought to be particularly relevant for magnetar QPO phenomena.) 
They found that the low current density approximation (LCDA) is valid for at least half of all toroidal oscillation modes 
and analyzed the energy losses due to plasma outflow caused by these modes for cases where the size of the polar cap (the 
region on the stellar surface that is crossed by open magnetic field lines) is small, finding that the energy losses are strongly 
affected by the magnetospheric plasma. For oscillation amplitudes larger than a certain critical value, they found that energy 
losses due to plasma outflow were larger than those due to the emiss ion of the electromagnetic waves (assuming in that case 
that the star was surrounded by vacuum). Recently. ITimokhin! (|2007l ) considered spheroidal oscillations of a NS with a dipole 
m agnetic field, usin g the TBS formalism, and found that the LCDA again holds for at least half of these modes. Discussion 



Timokhin] 1 2007 ) also provided some useful insights into the role of rotation for the magnetospheric structure of oscillating 



NSs. 

The TBS model was a very important contribution and, to the best of our knowledge, remains the only model for the 
magnetosphere of oscillating NSs available in the literature. However, it should be pointed out that it does not include several 
ingredients that a fully consistent and realistic model ought to include. Most importantly, it does not treat the magnetospheric 
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currents in a fully consistent way: although it gives a consistent solution for around half of the oscillation modes, the remaining 
solutions turn out to be unphysical and, as TBS pointed out, this is a symptom of the LCDA failing there. Also, rotation 
and the effects of general relativity can be very relevant; in particular, several authors ha ve stressed that using a Newtonian 
approach may not give very good results for the stru cture of NS magnetospheres (see, e.g., Beskin 199Cjl : Muslimov fc Tsvean 
19921 : iMofiz fc Ahmedovlbood : iMorozova et al.ll2008h . However, a more realistic model would naturally be more complicated 



than the TBS one whose relative simplicity can be seen as a positive advantage when using it as the basis for further 
applications. 

The aim of the present paper is to give a general relativistic reworking of the TBS model so as to investigate the effects 
of the changes with respect to the Newtonian treatment. We derive the general relativistic Maxwell equations for arbitrary 
small-amplitude oscillations of a non-rotating spherical NS with a generic magnetic field configuration and show that they can 
be solved analytically within the LCDA as in Newtonian theory. We then apply this solution to the case of toroidal oscillations 
of a NS with a dipole magnetic field and find that the LCDA is again valid for at least half of all toroidal oscillation modes, 
as in Newtonian theory. Using an improved formula for the determination of the last closed field line, we calculate the energy 
losses resulting from these oscillations for all of the modes for which the size of the polar cap is small and discuss the influence 
of GR effects on the energy losses. 

The paper is organized as follows. In Section 2 we introduce some definitions and derive the quasi-stationary Maxwell 
equations in Schwarzschild spacetime as well as the boundary conditions for the electromagnetic fields at the stellar surface. 
In Section 3 we sketch our method for analytically solving the Maxwell equations for arbitrary NS oscillations with a generic 
magnetic field configuration. In Section 4 we apply our formalism to the case of purely toroidal oscillations of a NS with a 
dipole magnetic field and also discuss the validity of the LCDA and the role of GR effects. In Section 5 we calculate the energy 
losses due to plasma outflow caused by the toroidal oscillations. Some detailed technical calculations related to the discussion 
in the main part of the paper are presented in Appendices A-C. 

We use units for which c = 1, a space-like signature (—,+,+,-1-) and a spherical coordinate system (t,r,0,<j>). Greek 
indices are taken to run from to 3 while Latin indices run from 1 to 3 and we adopt the standard convention for summation 
over repeated indices. We indicate four-vectors with bold symbols (e.g. u) and three-vectors with an arrow (e.g. u). 



2 GENERAL FORMALISM 

2.1 Quasi-stationary Maxwell equations in Schwarzschild spacetime 

The study of electromagnetic processes related to stellar oscillations in the vicinity of NSs should, in principle, use the coupled 
system of Einstein-Maxwell equations. However, such an approach would be overly complicated for our study here, as it is for 
many other astrophysical problems. Here we simplify the problem by neglecting the contributions of the electromagnetic fields, 
the NS rotation and the NS oscillations to the spacetime metric and the structure of the NS^i noting that this is expected 
to be a good approximation for small-amplitude oscillations. Indeed, for a star with average mass-energy density p, mass M 
and radius R, the maximum fractional change in the spacetime metric produced by the magnetic field is typically of the same 
order as the ratio between the energy density in the surface magnetic field and average mass-energy density of the NS, i.e., 

B2 io- f-?-V fM*®v_R_y . ( 4) 

mis p \ »,f m 1,™ / v / 



87rpc 2 V 1015G / V M /V 10km , 

The corresponding fractional change in the metric due to rotation is of order 

fi \ 2 ,„_,/ 11 \ 2 /l kHz 



where SIk is the Keplerian angular velocity at the surface of the NS. Moreover, in the case of magnetars, which we consider in 
our study, the oscillations are thought to be triggered by the global reconfiguration of the magnetic field. Due to this reason, 
the corrections due to the oscillations should not exceed the contribution due to the magnetic field itself given by estimate 
Q. Therefore, we can safely work in the background spacetime of a static spherical star, whose line element in a spherical 
coordinate system (t, r, 0, <j>) is given by 

ds 2 = goo(r)dt 2 + gn(r)dr 2 + r 2 d0 2 + r 2 sin 2 6d<j) 2 , (6) 

while the geometry of the spacetime external to the star (i.e. for r R) is given by the Schwarzschild solution: 

ds 2 = -N 2 dt 2 + N~ 2 dr 2 + r 2 d6 2 + r 2 sin 2 6dcj> 2 , (7) 



Several authors have, in fact, studied the equilibrium configurations of magnet ars by solving the Einstein-Maxwell e quations in full 
general relativity dBocqiiet et a l. 1995; Bonazz ola et al.l fl996: Cardall ct al.ll200lf ) or by using perturbative techniques ijColaiuda et al.l 
120081: lHaskell et al.ll2Q{m 
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where N = (1 — 2M/r) 1 ' 2 and M is the total mass of the star. For the part of the spacetime inside the star, we represent the 
metric in terms of functions A and as 

2#(r) 2A(r) /, 2m(r)\ 1 
3oo = -e ffii = e k; = I 1 ^-M , (8) 

where m(r) = 4-zr J Q r r' 2 p(r')dr' is the volume integral of the total energy density p(r) over the spatial coordinates. The form of 



these functions is given by solution of the standard TOV equations for spherical relativistic stars (see, e.g.. lShapiro fc Teukolskv 
19831 ) and they are matched continuously to the external Schwarzschild spacetime through the relat 



tions 



ff oo(r = R) = Nl , gu(r = R) = N~ 2 , (9) 

where N u = (1 — 2 A/ / R) 1 ^ 2 . Within the external part of the spacetime, we select a family of static observers with four- velocity 
components given by 

(/LeF 1 1,0,0,0). (10) 



and associated orthonormal frames having tetrad four vectors {efi} = (eg, ef , eg, e i) and 1-forms {u^} = («°, oj t , uj , w^), 
which will become useful when determin ing the "physical" componen ts of the electromagnetic fields. The components of the 
vectors are given by equations (6)- (9) of lRezzolla fc Ahmedov (h ereafter Paper I). 



The general relativistic Maxwell equations have the following form IjLandau fc Lifs hitz 1987) 

3-F]ti/3,7] = -Fq/3,7 + F-ya,P + Fl3~f,a = , (11) 

F a ^ ,p — 4irJ a , (12) 

where F a/3 is the electromagnetic field tensor and J is the electric-charge 4-current. We consider the region close to the star 
(the near zone), at distances from the NS much smaller than the wavelength A = 2-kc/lu. In the near zone the electromagnetic 
fields are quasi-stationary, therefore we neglect the displacement current term in the Maxwell equations. Once expressed in 
terms of the physical components of the electric and magnetic fields, equations (|lip and (|12p become (see Section 2 of Paper 
I for details of the derivation) 



sin0d r (r 2 B f ) + N^rdg (sin(9B e J + N^rd^B* = , (13) 

(rsin<9)^ =N [d^-dg (sinflE*)] , (14) 

(7V _1 r sin 6) = -d^E f ' + sin 6d r (rNE^ , (15) 

(ATV) 2|_ = (rNE § ) + d B E f , (16) 



N 



sin 9d r [r 2 E f ^ + rde (sin 6E 9 ^ + rd^E^ = Anp e r 2 sin 9 , (17) 



[de (smOB^ - c^B 8 ] = 47rr sin df , (18) 



d<t,B r — sin 8d r 



(rNB^ = 47rrsin(9J 8 , (19) 

d r (a^B 8 ) - dgB f = lirrJ^ , (20) 
where p e is the proper charge density. We further assume that the force-free condition, 

E S c ■ B = , (21) 

is fulfilled everywhere in the magnetosphere, implying that the magnetosphere of the NS is populated with charged particles 
that cancel the longitudinal component of the electric field. The charge density p sc responsible for the electric field E sc (cf. 
equation 1 1 7|) is the characteristic charge density of the force- free magnetosphere; this is appropriate for describing the charge 
density in the inner parts of the NS magnetosphere. We will refer to -B sc as the space-charge (SC) electric field, while to p sc 
as the SC charge density. 

Finally, we introduce the perturbation of the NS crust in terms of its four-velocity, with the components being given by 
a _ -*(, dx l \ -i>( , -Aj f Sv e Sv^ \ 

w = e V'Ht =e 1,6 5v ~ >7sm) • (22) 
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where 8v l = dx 1 /dt is the relative oscillation three-velocity of the conducting stellar surface with respect to the unperturbed 
state of the star. 

2.2 Boundary conditions at the surface of star 

We now begin our study of the internal electromagnetic field induced by the stellar oscillations. We assume here that the 
material in the crust can be treated as a perfect conductor and the induced electric field then depends on the magnetic field 
and the pulsational velocity field according to the following relations (see Paper I for details of the derivation): 



E\ n = -e~ 9 \5v B' p - 8v*B^ , (23) 
EL = -e~* \5v*B* - Sv^B^} , (24) 



Ef n = -e~* [Sv^B 9 - 6v 6 B*\ . (25) 

Boundary conditions for the magnetic field at the stellar surface (r = R) can be obtained from the requirement of 

continuity for the radial component, while leaving the tangential components free to be discontinuous because of surface 
currents: 

B ex \ r= n — B in \ r=R , (26) 

Bt\ r =R = Bi\ r=a + 4^ , (27) 

Bi\ r=R = Bi\ r =R-4m § , (28) 

where € is the surface current density. Boundary conditions for the electric field at the stellar surface are obtained from 
requirement of continuity of the tangential components, leaving E r to have a discontinuity proportional to the surface charge 
density E s : 

Zf ex \ r=R = E? n \ r=R + 4ttE s = -iV^ 1 [Sv^E* - 5v*B § ^ \ r=R + 4ttE s , (29) 

Ei\ r=R = Ei\ r=R = -iV^ 1 [Sv^B' - 8v f B^ \ r=R , (30) 

Et\ r=R = E? n \ r=R = -Nn 1 [Sv'B § - 5v § B f '] \ r=R , (31) 
where E s is the surface charge density. 

2.3 The low current density approximation 

The low current density approximation was introduced by TBS, and in the present section we present a brief introduction to 
it for completeness. Close to the NS surface, the current flows along the magnetic field lines, and so in the inner parts of the 
magnetosphere it can be expressed as 

J = a(r,e,<t>)-B, (32) 

where a is a scalar function. The system of equations (|13[l - (|20p . (|2ip and (|32[1 forms a complete set but is overly complicated 
for solving in the general case. However, within the LCDA these equations can, as we show below, be solved analytically for 
arbitrary oscillations of a NS with a generic magnetic field configuration. 

The LCDA scheme is based on the assumption that the perturbation of the magnetic field induced by currents flowing 
in the NS interior is much larger than that due to currents in the magnetosphere, which are neglected to first order in the 
oscillation parameter £ = £/R: 

4-7T -* -* 

-J«VxB, (33) 

c 

and 

V x B (1) = , (34) 
where B^ is the first order term of the expansion in £. This also implies that the current density satisfies the condition 
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J<< K S(0) !) c ^ sc(i?)c (^) ' (35) 

where p sc (R) is the SC density near to the surface of the star. Here we have used the relation p 3C (R) ~ B^°'rj/cR, where r\ 
is the velocity amplitude of the oscillation and u> is its frequency. 

In regions of complete charge separation, the maximum current density is given by p sc c. Since the absolute value of p sc 
decreases with increasing r and because r <C c/w in the near zone, condition (|35[) is satisfied in the magnetosphere if there is 
complete charge separation there. Since the current in the magnetosphere flows along magnetic field lines, its magnitude does 
not change and so condition (|35[) is also satisfied along magnetic field lines in non-charge-separated regions as long as they 
have crossed regions with complete charge separation. 

In the following, we solve the Maxwell equations assuming that condition (|35|) is satisfied throughout the whole near zone. 
As discussed above, a regular solution of the system of equations (|13|1 - (|20[) . (|21[) and (|32l) should exist for arbitrary oscillations 
and arbitrary configurations of the NS magnetic field and so, as shown by TBS, if a solution has an unphysical behaviour, 
this would imply that the LCDA fails for this oscillation and that the accelerating electric field cannot be screened only by 
a stationary configuration of the charged-separated plasma. In some regions of the magnetosphere, the current density could 
be as high as 

For a more detailed discussion of the LCDA and its validity, we refer the reader to Sections 2.3 and 3.2.1 of TBS. 



3 THE LCDA SOLUTION 

3.1 The electromagnetic field in the magnetosphere 

We now begin our solution of the Maxwell equations, assuming that the LCDA condition (|35[) is satisfied everywhere in the 
magnetosphere. Within the LCDA, equations (|18|l - (|20| l for the magnetic field in the magnetosphere take the form 

do (sin 05^ - d^B e = , (37) 

d4,B f - sin6»a r (rNB^J = , (38) 

d r (a^B 9 ) - deB" = . (39) 

As demonstrated in Paper I, the components of the magnetic field B r ,B and B^ can be expressed in terms of a scalar 
function S in the following way: 

B * = T^TR [sin Ode (sin 6d e S) + d^S] , (40) 

r z sin 6 

B § = — OedrS , (41) 
r 

B^ = -^L-d^drS . (42) 
r sm a 

Substituting these expressions into the Maxwell equations (|f 4p - (|16|) . we obtain a system of equations for the electric field 
components which has the following general solution 

E r sc =-d r (y sc ) , (43) 
E L=-J^d t d,S-^ sc ), (44) 

^^W-^M**,), (45) 

where ^/ sc . is an arbitrary scalar function. The terms proportional to the gradient of 9 SC are responsible for the contribution 
of the charged particles in the magnetosphere. The vacuum part of the electric field is given by the derivatives of the scalar 
function S. Substituting (|43[) - (|45[) into equation (|17[) . we get an expression for the SC charge density in terms of ty^: 



Psc = 



Nd r {r 2 d r t>sc) + ^Arf SC 



47]- r 2 

where An is the angular part of the Laplacian: 



(46) 
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An = -^—:de (sm8d e ) + -A-^d^. (47) 
sin (9 sir 9 



3.2 The equation for * s 



Substituting expressions (|40|| - (|42|) and (|43|l - (|45[) for the components of the electric and magnetic fields into the force-free 
condition (|21[) . we get the following equation for ^ sc 



-At- [sinOde (sm9d e S) + d^S] d r (* sc ) - -A- 
sm^ 9 sm! 



dodrSde{* 



d^dtSdedrS - dedtSd^drS] 
1 



d^drSdj, (* sc ) = 



(48) 



If the amplitude of the NS oscillations is suitably small (£ <C 1), the function S can be series expanded in terms of the 
dimensionless perturbation parameter £ and can be approximated by the sum of the two lowest order terms 



S(t,r,9,<f>) = S (r,9,<j,) + 5S(t,r, 



(49) 



Here the first term So corresponds to the unperturbed static magnetic field of the NS, while SS is the first order correction 
to it. At this level of approximation, equation (|48[) for ^g,-. takes the form 

1 



sin^fl 



i9d e (sm8dgS ) + d&So] ^(* sc ) - — , [d 4 ,d t (8S)ded r S - d dt(SS)d 4 ,d r S o 



dedrSode(V sc ) - — 
sin 



-fl^ft.5 e#(* BO ) = 



Next we expand S in terms of the spherical harmonics: 

CO I 

S = Y,Y. S lm {t,r)Y lm {9,4>) . 



(50) 



(51) 



where the functions St m are given in terms of Legendre functions of the second kind Qi by (jRezzolla et al.ll200lf ) 



. . r d 



Stm{t) 



(52) 



Note that all of the time dependence in (|52|l is contained in the integration constants si m (t) which, as we will see later, are 
determined by the boundary conditions at the surface of the star. We now series expand the coefficients St m {t, r) and st m {i) 
in terms of £ 



Stm(t, r) = S im(r) + SSi m (t,r), s lm (t) = s olm + 5s tm (t) , 



(53) 



where all of the time dependence is now confined within the coefficients SStm(r,t) and 8se m (t), while the coefficients Soe m 
and soi m are responsible for the unperturbed static magnetic field of the star. Using these results, we can also express S and 
8S in terms of a series in Yi m (9, <f>) in the following way 



So = ^2 Soern(r)Ye m (9, 



SS = Y.Y. SS em (t,r)Y em (9, 



(54) 
(55) 



The variables r and t in the functions Se m (t,r) and Se m (t,r) can be separated using relation (f52 



Sotm{r) 



r 1 d 
M 2 dr 



r 2 d 



r J dr \ 

1-^)^(1 
r I dr \ 



r s 
~M / 

r 

Jl 



Ssim(t) 



(56) 
(57) 



3.3 The boundary condition for \l/ gc 

We now derive a boundary condition for ^ sc at the stellar surface using the behaviour of the electric and magnetic fields in 
that region. Following TBS, we assume that near to the stellar surface the interior magnetic field has the same behaviour as 
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the exterior one: 

B " = — T^T-n\ An6d <> [s^OdeS) + d^S\ , (58) 
r 2 sin 9 

P -A 

B s = Ci decks' , (59) 

r 



-A 

= Ci — : — -dcpdrS . (60) 



U sing the continuity condition for the normal component of the magnetic field [B r ] — at the stellar surface ( Pons fc Geppert 
2007 ) together with the condition e _A | r= ij = Nr, one finds that the integration constant C\ is equal to one. The interior 



electric field components can then be obtained by substituting ()58p _ <160p (with C\ = f) into (|23p - (125p : 

p -(*+A) , . i -, 

El n = tt- I 5v e dsd r S - sm96v rt, d e d r S \ , (61) 

rsinff I J 

= { Sv'd^drS + [sin 6de (sin 0<9 fl S) + } , (62) 

S » = -*— r {^vd e d r S + [sm9d e (smOdeS) + 8^S\ \ . (63) 

The continuity condition for the 9 component of the electric field across the stellar surface (|30|) gives a boundary condition 

scl 

Sv* 
Rsin 2 9 



for de^f sc \ r =R 

* sc 9 | r=H = - <{ [sin<9<9 e (sin 98 e S) + d^S] + i^-d^A-S + ^dtd^S \ \ r =R , (64) 



while the continuity condition for (|31 1) gives a boundary condition for sc \ t —r: 



*sc,Jf=H= i [ sm (sin 9d g S) + d^S] + NSv r sin 9d g d r S + sin 0c> t <9 9 S }> | r=jR . (65) 

xx Sill [/ 



Integration of equation (|64[) or equation (|65p over 9 or (f> respectively, gives a boundary condition for * sc . We will use the 
result of integrating equation (|64p over 9. Assuming that the perturbation depends on time t as e" 1 ^*, we obtain the following 
condition, correct to first order in £, 

( \ ^2 n l sindd <> (sin^So) + d4.d4.S0] + ^-d4,drS + J—d t d4,(5S) \ d9\ r=R + e iwt F((j>) , (66) 
J I H sin 8 sin 9 sin 9 

where F(cj>) is a function only of cj> which we will determine below. 

The components of the stellar-oscillation velocity field are continuously differentiable functions of r, 9 and <f>. The boundary 
conditions for the electric field (|30|l - (|3ip imply that the tangential components of the electric field E sc must be finite. The 
vacuum terms on the right-hand side of (|44p - (|45p and the terms on both sides of equation (|45[1 are also finite. Consequently, 
the term 

" -^9~ lr = R (67) 
should also be finite. Hence we obtain that d4,(<li sc )\g =0 n . r=R = and so the function F(<j>) in the expression for boundary 
condition (|66p must satisfy the condition ($ gc )|e =0j7r;r . = i{ = C e -1 ^*, where C is a constant. Using gauge invariance, we choose 

# sc |0=O;r=fl = O, (68) 

and from this and equation (|66[) . we obtain our expression for the boundary condition for ^/ sc at the stellar surface: 



x v <t> N/)v f 1 



4 TOROIDAL OSCILLATIONS OF A NS WITH A DIPOLE MAGNETIC FIELD 

As an important application of this formalism, we now consider small-amplitude toroidal oscillations of a NS with a dipole 
magnetic field. For toroidal oscillations in the (£',m') mode, a generic conducting fluid element is displaced from its initial 
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location (r, 9, </>) to a perturbed location (r, 9 + £ e , <j> + f^) with the velocity field ( Unno et al, 19891 ). 

Sv* = , fa* = 4- = e'^rjir) J-^F fw (9, 0) , 6v $ = = -e^Md^w (A, 0) , (70) 



where uj is the oscillation frequency and rj(r) is the transverse velocity amplitude. Note that in the above expressions (|70l> . 
the oscillation mode axis is directed along the z-axis. We use a prime to denote the spherical harmonic indices in the case of 
the oscillation modes. 



4.1 The unperturbed exterior dipole magnetic field 

If the static unperturbed magnetic field of the NS is of a dipole type, then the coefficients sof m involved in specifying it have 
the following form (see eq. 117 of Paper I) 



SO 10 = — - 



/3tt 



-/i cos x 



son = \j y/ismx 



(71) 



where /j, is the magnetic dipole moment of the star, as measured by a distant observer, and x is the inclination angle between 
the dipole moment and z-axis. Substituting expressions (|7ip into (|56p and then the latter into ()54p , we get 



S = 



3^r 2 
8M 3 



IniV H 



1 + ^ 
r 



(cos 9 cos x + e sin sin x) 



The corresponding magnetic field components have the form 



B r 



3fj, 
~4M 3 

AM 2 r 
3fj,N 



In AT H 1 



M 
;■ 



^-lnAr 2 + -L + l 



(cos x cos + sin x sin 0e ) , 
(cos x sin — sin x cos 0e ) , 
(-isinxe 10 ) . 



(72) 

(73) 
(74) 
(75) 



4MV 

At the stellar surface, these expressions for the unperturbed magnetic field components become 

B r R — /h-^o(cosxcos0 + sinxsinfe"*) , B R = h R B (cos x sin 9 — sin x cos #e"^) , .BJ^ = — i/i^B^sinxe 1 ^) , (76) 

where Bo is defined as Bo = 2/i/R 3 . In Newtonian theory Bo would be the value of the magnetic strength at the magnetic 
pole but this becomes modified in GR. The GR modifications are contained within the parameters 



h„ = 



3R 2 N B 
8M 2 



fn 



3R 3 
8M 3 



In Nz 



2M 



, M 
R \ 1+ R 



(77) 



Fo r a given ji, the magne t ic fiel d near to the surface of the NS is stronger in GR than in Newtonian theory, as already noted 
bv lGinzburg fc Ozernov I 1 19641 ). 



4.2 The equation for ^ sc 

Substituting So from (|72p into equation (|50p . we obtain a partial differential equation containing two unknown functions ^ sc 
and SS for arbitrary oscillations of a NS with a dipole magnetic field 

— 2r 2 ?i(r) (cos^cosx + e !< ^ sin^sinx^ 9 r (* sc ) + d r [r 2 qi(r)] (sin 9 cos x — ^ cos 9 sin x) ^e(^ sc ) (78) 



-d r [r 2 q!(r)] 



e' v sm x 



0*(* S c) + 



9. [r 2 qi(r)} 



[(sinScosx — e"^ cos 9 sin xj d^d t (SS) + ie"^ sin x sin #<9 9 <9 f (5S)J = 



sm # sm f 

where we have introduced a new function q\ (r) for simplicity of notation [see Eq. (|A2p for the definition of qi (r) 
From (|69p . the boundary condition for ^ sc at the stellar surface is 



SC \r = R 



B(,Rf R 8v^de (cos 9 cos x + e 1 ^ sin sinx 



sin 9 



d t d^8S)\d9\ r=R . 



(79) 



Using the expressions for the velocity field of the toroidal oscillations (|70p and for the boundary conditions for the partial 
derivatives of the SC potential (|64p ~ (|65p . we find that d t SS is given by (see Appendix A for details of the derivation) 



dtSS(r,t)=J2 



£(£+1) R 2 qe(R) 



(80) 
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I J 



dgYprn ( sin 9 cos x — ^ cos sin \ ) + ie"^d$Yt m sin 9 sin x\ 



y;, 

sin 9 



From here on, for simplicity, we will consider only the case with % — 0- Although our solution depends on the angle between 
the magnetic field axis and the oscillation mode axis, focusing on the case \ — does not actually imply a loss of generality 
because any mode with its axis not aligned with a given direction can be represented as a sum of modes with axes along this 
direction. We have developed a MATHEMATICA code for analytically solving equation (|50|) and hence obtaining analytic 
expressions for the electric and magnetic fields and for the SC density. 

The solution of equation (|T8[) for the case x = is given in Appendix [Bj where we show that the general solution has the 
following form 

,t, 1 ™' 2 Df _ f r d r ,[r' 2 qi (r')} g e ,(r') W(g(rQ^) , r > s - - . , ,1 

where r' is the integration variable. In order to solve this integral, the function 8(r') is expressed in terms r' and a constant <fi2 
through the characteristic equation (|B6[) and, after performing the integration, tfi2 is removed again using (|B6fl . The unknown 
function $2 is determined using the boundary condition for $ 2 |r=_R given by (|F313fl . Once the integral on the right-hand side 
of (EH has been evaluated, we then express all of the trigonometric functions resulting from the integral, in terms of sin 9. 
This 9 is the value at r — R. To obtain an expression for the value of $2 at a general radius, we write this 9 (at r = R) in 
terms of the value of 9 at a general point (with r > R) using the characteristic relation (|B6|l . i.e., 



sintf^ V[r 2 9iM]/[#V(#)] x sinS , (82) 
so that 

2 1 \ -1 1 / 2 

- fl2 gi (fl) sm 9 sign(cos 6>) , (83) 

where sign(a;) is defined such that sign(a;) = +1 if x > 0, and sign(a;) = —1 if x < 0. There are then different expressions 
for ^ sc in the two regions 9 € [0, 7r/2] and # G [tt/2, 7t]. If these two expressions do not coincide at the equatorial plane for 
r > R, then there will be a discontinuity in 9 SC at 9 = 7r/2, and quantities that depend on dgty sc will become singular there. 
As shown by TBS, the function 9 SC is indeed discontinuous at 9 = n/2 for some oscillation modes and, as we discussed in 
Section T2.3I above, this unphysical behaviour indicates that the LCDA ceases to be valid for those modes. In these cases, the 
accelerating electric field cannot be canceled without presence of strong currents which may become as large as (|36[l in some 
regions of the magnetosphere. The occurrence of such singularities was explained by TBS and the reader is referred to Section 
3.2 of their paper for a detailed discussion. 

Next we discuss how GR effects contribute to our solution. As discussed above, for a given magnetic moment fj, (as 
measured by a distant observer) the strength of the unperturbed magnetic field near to the surface of the NS is larger in GR 
than in Newtonian theory. Due to the linearity of the Maxwell equations, a perturbation of a stronger magnetic field should 
produce a larger electric field for the same oscillation parameters. This in turn should lead to a larger absolute value of the 
SC density in GR, since the SC density takes the value necessary to cancel the electric field. In the next Section, we will give 
a more quantitative analysis of the GR contribution in our solution. 

We point out that the function ^ sc does not depend on £ > 1 perturbations to the stellar magnetic field in the case 
of axisymmetric (m' = 0) toroidal modes. This can be seen from the fact that these perturbations are confined within the 
5S terms which enter equation (I78|l for the function ^ sc only through a derivative with respect to <f>; hence vanish for the 
axisymmetric modes. Therefore, the only perturbation to the magnetic field is due to the 1=1 term, and the solution for 
these modes is much simpler than that for non-axisymmetric (m' 7^ 0) modes. It is then convenient to discuss separately the 
axisymmetric and non-axisymmetric cases. 

The solution (|81|l for case the m! = modes at r — R has the following form 

* sc (r,M,t)|r= H = -B Rf R fj R f cGB0a*Yi» o (0,0 d-d . (84) 

Jo 

Using the properties of the spherical harmonics, we can express *t> sc (r,9,<f>,t) \ r =R for odd £' modes in the general form 

N 

* sc (r, 9, </>, t) \ r=a ~J2 A 2n sin 2 ™ 6 , (85) 

n - 1 

while for even £' modes, the general form of 9 SC is 

N 

* sc {r,9,cj>,t)\ T=R ~ (yl + Bcos60^yl2nSin 2n <9 , (86) 

n=l 

where the coefficients A, B and A2 n do not depend on r and 9. The value of N equals £' /2 + 1 for even £' and (£' + l)/2 for 
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r/R r/R 

Figure 1. Left panel: The ratio (^ sc )gr/(^sc ) Novrt along the polar axis plotted as a function of the distance from the star, for 
axisymmetric toroidal modes of a NS with compactness M/R = 0.2. Right panel: The ratio (p sc )gr/(Psc )Newt m * ne equatorial plane 
plotted as a function of r, for a star with M/R = 0.2, for toroidal oscillation modes (2, 0), (1, 1) and (3, 0). 



odd I'. As we discussed above, in order to obtain the solution for 9 SC for r > R, one has use sin# — > \J r 2 qi (r) / R 2 qi (R) sin# 
on the right-hand sides of (|85|) and (|86p . Thus the GR effects contribute to the solution for m! = modes only through terms 
fn [ T ' 2 qi ( r ) /R 2 Qi(R)] " > where n ^ 1. Note that the factor f R in this term appears due to the boundary condition at the surface 
of the star, namely from the continuity of the tangential components of the electric field, while the factor [r 2 qi(r)/R 2 qi (R)] ™ 
appears due to the presence of charged particles in the magnet osphere. The second factor is equal to 1 at the stellar surface 
and approaches its Newtonian value (R/r)" at large r and small M/R. Since f R > 1, the absolute value of 5' sc should be 
greater in GR than in Newtonian theory. For example, in the case of small 8, the only term which is important is that with 
n = 1 and hence we get (^ sc ) GR /(^' SC ) Nowt = f R r ' l <li( r )/R' i< li(R)- This quantity is shown in Figure Q] (left panel), where 
we can see that, near to the stellar surface, the function \I/ gc is larger in GR than in Newtonian theory, while at larger r it 
asymptotically approaches its Newtonian value. 

Analysis of the GR contribution to the solution in the case of non-axisymmetric modes is more complicated because in 
this case the solution depends not only on £ = 1 perturbations to the magnetic field but also on I > 1 perturbations, which 
are contained in the term dtd^SS of equation (I78|) . and contribute to the solution due to the integral in (|8ip . Nevertheless, 
some rough estimates of the GR effects can be made in the following way. Near to the stellar surface the integral in (|81[1 can 
be approximated as 

r d4z 2 qi (z)] g e (z) Y tm (0(z),<j>) r 2 g t {r) Y lm {6{r),d>) 

Ja J R gi (z) R 2 qi{R) cos6»0) Jr R 2 q e (R) cos(9(r) ' { ' 

Close to the star, (r 2 qi(r))/(R 2 qt(R)) ~ 1 and so the leading GR contribution comes from the factor f R which increases the 
absolute value of ^ sc with respect to the Newtonian case. Further away from the star, r 2 qi is approximately proportional to 
r~ e and so the integral in (|87[1 can be approximated as ~ .f R {R/ r Y +1+m ' 2 to leading order in R/r. Therefore, while this 
integral makes an important contribution to ^ sc near to the star, it becomes negligibly small for £ > 1 at r > R as compared 
with $2- The GR effects contribute to $2 through the terms f R [r 2 qi(r) / R 2 qi(R)^" in a similar way to their contribution to 
^ sc for the axisymmetric modes discussed above. This increase of the function ^ sc due to the GR effects also lead to an 
increase in the absolute values of the SC density p sc near to the star, as shown in Figures [T] (right panel) and [2] for some 
toroidal oscillation modes. 



5 ENERGY LOSSES 

It was shown by TBS that the kinetic energy of the stellar oscillations should be lost through being passed to plasma near 
to the stellar surface which then flows out along the open magnetic field lines. Note that within the framework of the TBS 
model, electromagnetic fields are considered only in the near zone and so the existence of the plasma outflow cannot be shown 
explicitly; however, qualitatively, the mechanism for the plasma outflow should be the following. The charged particles that 
were accelerated to high energies by the longitudinal electric field move along the magnetic field lines in the near zone. If the 
kinetic energy density of the plasma at the equator becomes comparable to the energy density of the magnetic field at some 
point, then the field line which crosses the equator at that point becomes open. Plasma flowing along open field lines forms 
an electromagnetically driven wind which closes at infinity. There is then an electric current flowing along the stellar surface 
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Figure 2. The space-charge density p sc at r = 1.5R (in units with Bq = rj = 1) plotted as a function of 9 for the toroidal oscillation 
modes (1,0), (1, 1) and (3,0). The solid lines refer to the SC density p sc in the relativistic theory, while its Newtonian value is shown 
by the dashed lines. 



between positive and negative emission regions. Because this current must cross the magnetic field lines at the stellar surface, 
it exerts a braking torque on the NS oscillations and thus reduces their kinetic energy (see Section 3.2.2 of TBS for more 
details). 

In the following, we carry out a GR calculation of the energy lost by toroidal stellar oscillations due to plasma outflow. 
First we calculate the energy losses due to the outflow of a particle along an open field line from a given point on the stellar 
surface. For this purpose, we start by considering the motion of the charged particle along the ^-direction on the stellar surface 
(where it crosses the magnetic field lines, hence exerts a braking torque on the stell ar oscillations). The equ ation of motion 
for a test particle of mass m in a generic electromagnetic field has the general form ( Landau fc Lifshit3 1987 ) 



m — ; — = eb H WQ , 
dr 

where D/dr is a comoving derivative, w " is the four- velocity of the particle given by 
u a + v a 



VT 



(88) 



(89) 



u a is the 4- velocity of the static observer (|10[) . and d° is the velocity of the particle relative to the static observer. 

Because of time-invariance, there exists a timelike Killing vector £° such that = —N 2 . The four- velocity of the static 
observer can be expressed in terms of £ a as u a = N^ 1 ^", and therefore the energy of the particle is given by 



I = — p a £ = —mvjat, . 
Contracting the equation of motion 

ZamWtfW^ = -eF M £a , 

Vl — v- 1 

The right-hand side of this can be rewritten as 



with the Killing vector £ a gives 



eb £fj — eb H - 



= eF ap N- Uf3V ° 



E a v&N 



while the left-hand side can be transformed as 

j- a 3 I j-ol \ 3 ct 3 

^ciriw .pw — m(w a t; );3W — mt, a -fsw w , 



(90) 



(91) 



(92) 



(93) 



The second term on the right-hand of this equation vanishes due to antisymmetry of the tensor f a ;/3- Therefore, the projection 
of the equation of motion onto the Killing vector can be written as 



dS 

717 



eN 



E a v & 



VT 



For particles moving along the ^-direction, this equation takes the form 

E & RN 2 R JO 
VI -v 2 



(94) 



(95) 
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Integrating this over the interval (0,0), we get 

AS = eRN R / E a d6. (96) 



Note that in deriving this last equation, we have used the fact that v <C 1. The quantity AS measures the energy that would 
be carried away by a particle that leaves the surface of the star from a point with coordinates (R, 9, <f)). In order to calculate 
the energy loss per unit time through a given surface element dS, AS needs to be multiplied by the the current density, j a 
and integrated over the surface element n a dS, where n a is the unit spacelike vector orthogonal to the surface r — R. The 
energy loss per unit time due to plasma emission through the surface element dS is then given by 

dL = AS(6,(f>)j & (R,0,(t>)na R 2 sin 6d8dtj> , (97) 

while the total energy loss L from the NS is obtained by integrating dL over the entire open field line region on the stellar 
surface. 

Now we determine the angle 8o at which the last closed magnetic field line intersects the stellar surface. Following TBS, 
we define the last closed line as being that for which the kinetic energy density of the outflowing plasma at the equator 
becomes equal to the corresponding energy density of the NS ma gnetic field. We now der ive a mathematical condition for 



this. The energy-momentum tensor of the electromagnetic field is ( Landau fc Lifshita 1987 ) 



T:l = ± [F^F? - Ig^F^F^j . (98) 
Using this expression, one can obtain an expression for the energy density of the electromagnetic field 

S em = NT al3 u a iip = N B2 + E2 . (99) 

The method for the calculation of the last closed field line proposed by TBS [equation (65) of TBS] is based on an 
implicit assumption that the plasma flows out from the star isotropically and its kinetic energy is distributed uniformly over 
the surface of a sphere with radius r = R a , where R a is the radial coordinate of the point where the last closed field line 
crosses the equatorial plane. However, according to the definition of the last closed field line, the outflowing plasma should 
move along the field lines throughout the region r < R a and so its energy cannot be distributed uniformly over the sphere 
r = R a . In the following, we derive an alternative formula for the calculation of do that takes into account the correction due 
to the anisotropic plasma outflow along the magnetic field lines. 

If the stellar magnetic field is dipolar, then the density of the field lines decreases with distance from the centre of the 
star. Since the charged particles move along the field lines, this implies that the kinetic energy density of a comoving element 
of plasma should decrease monotonically with r as the field lines separate (we are assuming here that the plasma speed along 
the field lines is roughly constant). Since for dipole magnetic field lines f r sva?6/r = const, it is easy to show the energy 
density of the plasma when it reaches the equatorial plane is smaller by a factor of 0.5R 3 /(R a f R ) compared to its value at 
the stellar surface (assuming that the angle 8 is small). Therefore, the kinetic energy density of the outflowing plasma at the 
point (R a , 0, <j>) is given by 

e pl (R a ,0,4>) = ±.£jL?L.A£j*n tt , (100) 



where Nu a = \/l — 2M/R a and the factor Nr/Nr^ accounts for the gravitational redshift of the energy of the plasma. The 
energy density of the magnetic field lines at the same point is determined from (|99|l and has the following form 



Nr R i 

te m (Ra,O,0) = -^—rB 2 o , (101) 
OZ7T H a 

Note that R a ^> R for small 80, and so in this case N R w 1 to good accuracy; from here on we will consider only the case of 
small 6*o and take N R = 1. Also j r 3> j 6 for small 8, and so in the following we will neglect the 8 component of the current 
above the polar cap. The last closed field line is determined by the conditions 

e p i(Ra,0,(f>) = e em (R a ,0,(f>) ■ (102) 
Substituting (|101|) and (|100[) into (| 102|) . one can obtain an algebraic equation 

16^iVflA£:p sc = f 4 R B 2 8t (103) 

for determining 80 after expressing AS in terms of 80 . Once we know 80, we can calculate the total energy of the outflowing 
plasma by integrating (|97|) over the entire open field line region of the stellar surface: 

d8 \f(R,9,(j>) AS (9, <j>)\R 2 sin8 . (104) 
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The motion of charged particles in strong magnetic fields can be approximated as a relativistic motion along the field lines. 
This is a reasonable approximation because the electrons are estimated to beco me relativistic at a height o f a few centimeters 
above the stellar surface. Therefore in the frame of the static observer we take (|Muslimov fc Tsygarj[l992T) 

/ = pc^ (105) 

where p(R, 8, (f>) is the SC density. 

TBS solved the Newtonian version of (|103[F1 and calculated L for three modes: (1,1), (2,0) and (3,0). They also estimated 
the order of magnitude of the energy loss L and the angle 8o for all of the oscillation modes for which 8o is small. They found 
that 8o is small for modes with m! < 3 and, as we will see later, this result also holds in GR. We will next solve the relativistic 
equation (|103p and calculate 8o and L for all of the modes (£\m') for which the angle 8o is small. We will present elsewhere 
a study of the case for large (9o. 

In order to solve equation (|103p and calculate L in the case of small 8q, we use the following procedure. First, we 
approximate equation (|B2|I and the boundary condition (|B1|) for small 8 by expanding all of the functional dependence on 8 
in a Taylor series, taking into account only terms of the two lowest orders in 8. Secondly, we solve this equation and calculate 
p sc and A£ using the approximate solution for ^f sc . Then, we calculate the angle 8q by substituting p sc and A£ into (|103[) . 
Finally, substituting the expressions for do, p sc and A£ into (|104[) . gives us the expression for L. Since the methods for solving 
equation (|B2|) for m! = and m! 7^ exploit a similar technique, we present just the solution for m' = below, in this Section, 
while the solution for non-axisymmetric modes is presented in Appendix [C] 

For small angles 8, the spherical harmonic Yi m {8, <j>) can be approximated by the sum of the first two lowest order terms 
in the expansion in terms of 8 

Y lm {8, <f>) « 8 m + A<*{4>) e m+2 ■ (106) 

In this case, the characteristics (IB6I) take the form 



tp= ^J-r 2 qi {r)8 . (107) 
Specialising to small 8 in the boundary condition (|B1[1 for $ gc , gives the following expression: 

9 aa \r=* = -B Q Rf R f, R Af 9 2 . (108) 



In order to calculate V&sc for arbitrary r, one has to replace 8 in this equation by \/[r 2 q\ (r)]/[R 2 q\ (R)] 8, following the 
argument given earlier for sin 8, and this then gives the solution: 

* sc = -B Rf R fj R jfi^Af Q 8 2 . (109) 



Substituting this into formula (|46(l . we obtain the expression for the SC density: 



BofnVR r 2 qi(r) a{2) 
n N R R R 2 qi(R)' 



psc = „zr^ ■ (no) 



Substituting ^sc as given by (|109|l into (|44|l . we obtain the expression for electric field E^„. Substituting this into (|96[) gives 
AS: 

A£(8, 4>) = B Rf R ^ R N Ti A (2 l 8 2 . (Ill) 
Using this expression for A£ in equation (|103|l . we obtain an algebraic equation for #o, which has the following solution: 



8 = 2NI/ 4 



r> A (2) 



1/2 



(112) 



fn 

Substituting (|llip and (|105p into (|104p . gives us the following expression for the total energy loss: 

Le = \ [B Rf R ii R A$l\ 2 8 4 , (113) 
and then using (|112l) in equation (|113l) . gives 

L eo = 8N R \b R f R (A {2) ) 2 ] * . (114) 



2 It should be noted here that, as we mentioned above, TBS uses a slightly different formula for the determination of 9q which does not 
take into account the anisotropy of the plasma outflow. 
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The energy losses for non-axisymmetric modes are calculated in Appendix [C] The angle Go and the energy loss L are given by 
expressions ()C7[1 and ()C9|) in that case. Note that determining 80 using the direct relativistic extension of the TBS method 
would give a value that is smaller than (|112[) by a constant factor of 2 ' 5 , causing L to be smaller by a factor of 4. For the 
ml = 1 modes, these correction factors are (27r) ' 25 and 2-zr respectively, while for the ml = 2 modes they are (37r) ' 5 and 

(37T) 3 . 

We now continue our discussion with analysis of the contribution of the effects of GR to the quantities 6q and L e i . The 
ratio of 9 in GR to its Newtonian counterpart for a given oscillation amplitude [i.e., for {fj R ) GR = (^ K ) Nowt ] can De obtained 
using equation (|112[) : 

(do) N 1/4 

v 'Newt JR 

This ratio is shown in Figure (left panel) plotted as a function of the stellar compactness parameter M/R. It can be seen 
that in GR, the angle 60 is smaller than in Newtonian theory (it approaches its Newtonian value as M/R — > 0), meaning 
that GR effects lead to the polar cap being smaller. The reason for this can be understood in the following way. The angle 
60 is determined from equation (|103[) and it is convenient to analyze the relativistic effects by looking at this equation. The 
contribution from the curvature is contained here in three terms: i) the factor on the right-hand side, u) the factor Nr 
on the left-hand side and Hi) a factor on the left-hand side that is contained implicitly within the term A£p sc . The first 
of these three terms accounts for the modification of the geometry of the background dipole magnetic field lines due to the 
curvature. As we discussed above, the magnetic field lines obey the relation f r sinO/r 2 = const, with f r = 1 in the Newtonian 
case but f r > 1 in GR. This means that a magnetic field line with a given energy density at r = R a on the equatorial plane, 
crosses the stellar surface closer to the pole in GR than in Newtonian theory. The second relativistic factor accounts for 
the redshift of the kinetic energy of the plasma due to the displacement from r — R to r — R a . Finally, the third factor is 
responsible for the amplification of the energy density of the outflowing plasma caused by the increase of the magnetic field 
strength (for a given magnetic moment) due to GR effects. The last of these would obviously lead to an increase of #0 in GR 
if it were acting alone. However, this effect is counteracted by the change of the geometry of the magnetic field lines and the 
gravitational redshift and those two effects are substantially stronger, giving the result that the angle 9o is smaller in GR than 
in Newtonian theory. Note that, as expected, the ratio (#o) GR / (^o) N for the axisymmetric toroidal modes does not depend 
on £' , while for the general ml 7^ modes, there is dependence on £' . This is because, as discussed above, the function ^ sc 
for m' — depends only on the lowest £ = 1 perturbation to the magnetic field since the higher I perturbations are contained 
in the term d^dtSS which disappears in axisymmetry. 

The ratio of the angle 60 in GR to its Newtonian counterpart, for ml 7^ modes, can be obtained using equation (|C7|l 
(see Appendix [Cl for details of the derivation): 



( fl oL B (N R \D« 



'I'm' Igr 



(^) N8wt V Pr l-Df'm'Lcwt 



(116) 



Figure [3] (central and right panels) shows the ratio of angle 60 in GR to its Newtonian equivalent for the modes ml — 1,2 
for several values of £'. As can be seen in these plots, 9q is smaller in GR than in Newtonian theory also for the ml = 1 and 
ml — 2 modes. The reason for this is similar to that for the ml = modes which we discussed above. We have made the 
analysis for values of £' up to t' = 10 and the ratio (#o) GR / (#o) N was found to be rather insensitive to the values of £' . 
(The dependence on I 1 is due to the fact that (p sc ) GR / (p S c) N t ^ s generally larger for higher £' at small 8o-) 

The ratio of the energy losses for the m' = modes in GR and in Newtonian theory is equal to Nr, as one can obtain 
using equation (|114j) . and this is plotted in the left panel of Figure [4] The GR modification is caused by the gravitational 
redshift of the plasma energy; the other "magnetic" GR effects do not influence this, because the shrinking of the polar cap 
and the increase in the plasma energy density exactly compensate each other. 

The ratio of the energy losses for the m' 7^ modes in GR and in Newtonian theory can be straightforwardly obtained 
using equation (|C9|l : 

=f R Np f/- 2 ,^; 1 - )^ . (nr) 

l-k<"m'J Nowt V \ L> l'm' I Newt / 

This is also shown in Figure [4] where one can see that the energy loss of the ml 7^ modes is smaller in GR than in Newtonian 
theory. The reason for this is the same as for the ml — modes discussed above, i.e., the increase in the energy density of 
the outflowing plasma cannot compensate the shrinking of the polar cap. Moreover, for the ml = 2 modes, the energy density 
of the outflowing plasma is proportional to 8 2m (see eqs. (|96[> and l|C5|l ) and so the total energy losses are very sensitive to 
the size of the polar cap. Because of this, the ratios of the total energy losses for the ra' = 2 modes in GR and in Newtonian 
theory are much smaller than those for the ml = and ml = 1 modes. 
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Figure 3. The ratio (9q) gr /(^o) New t * s plotted against stellar compactness M/R for some representative toroidal oscillation modes with 
m! = 0, 1, 2; it asymptotically approaches 1, of course, in the Newtonian limit M/R — » 0. 
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Figure 4. The ratio L Gn / L NBvrt is plotted against stellar compactness M/R for toroidal oscillation modes with m' = 0, 1, 2. 



6 SUMMARY 

In this paper, we have described our general relativistic model for the force-free magnetosphere of an oscillating, non-rotating 
neutron star. Our approach is based on the previous Newtonian model developed by TBS and focuses on toroidal modes 
which are thought to be particularly relevant for magnetar QPO phenomena. We have taken the spacetime geometry to be 
spherically symmetric and have neglected any modifications of it caused by the electromagnetic fields, the stellar oscillations 
and the magnetospheric plasma. Within this context, we have derived the general relativistic Maxwell equations for arbitrary 
small-amplitude oscillations of a neutron star with a generic magnetic field configuration and have shown that, as in the 
Newtonian case, they can be solved analytically for the force-free configuration of the electromagnetic fields (i.e. <C B ) 
under the low current density approximation (LCDA) . We have applied our formalism to small-amplitude toroidal oscillations 
of a neutron star with a dipole magnetic field and have found that the LCDA is valid for at least half of these modes in GR, 
as in the Newtonian calculations of TBS. We have also discussed the contribution of the GR effects to our solution, finding 
that they lead to an increase in the absolute values of the electromagnetic fields and the space charge (SC) density near to 
the stellar surface. 

We have calculated the energy losses due to plasma outflow resulting from these oscillations, focusing on cases where the 
size of the polar cap is small so that one can expand the Maxwell equations as Taylor series in powers of 8, retaining only the 
two lowest order terms. This approach leads to a great simplification and allowed us to perform a thorough analysis of the 
solution. We have found that in GR, the polar cap is smaller than in Newtonian theory and have shown that this is due to 
the change in the geometry of the dipole magnetic field and the gravitational redshift of the energy of the outflowing plasma. 
Also, we found that the oscillation modes which have small 8o, all have m' < 3 as in the Newtonian case. 

The total energy loss resulting from the stellar oscillations causing plasma outflow through the polar cap region, is 
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determined through an integral over the whole polar cap area, and so it depends on both the kinetic energy density of the 
outflowing plasma and the surface area of the polar cap. Although GR effects lead to some increase in the energy density of 
the outflowing plasma (due to the increase in the surface magnetic field strength for a given magnetic moment), the area of 
the polar cap is smaller in GR and we have found that the increase in the energy density of the outflowing plasma cannot 
compensate for the shrinking in size of the polar cap. Therefore the total energy losses for the toroidal oscillation modes are 
significantly smaller in GR than in Newtonian theory. 

In conclusion, we point out that while our calculations represent an advance with respect to previous ones, they still do 
not include a number of very important aspects which would be necessary for a realistic description of these phenomena. Most 
importantly, as noted above, they do not take account of electric currents flowing in the magneto sphere. Inclusion of these in a 
consistent way would require solving a version of the non-linea r "pulsar equation" ( Michel 1 19731 ) fo r oscillating neutron stars 



possi bly by adopting a numerical approach similar to that of IContopoulos et al. ( 199S ) (see also 



_Ljl_ 

20061 ) or performing time- dependent simulations of the magnetosphere ( Spitkovskvl 2006 ; Komissarov 



Gruzinov 
2006|: 



20051; iTimokhin 
McKinnevlbooel ) 



This will be the subject of future investigations. 



ACKNOWLEDGMENTS 

We wish to thank Luciano Rezzolla for bringing this problem to our attention and for many helpful comments. We also 
acknowledge helpful discussions with V. S. Beskin, U. Geppert, K. Glampedakis, A.V. Khugaev and A.N. Timokhin. BJA 
gratefully acknowledges travel support from the Exchange of Astronomers Programme of the IAU C46-PG-EA and the ICTP 
Associate Programme. 



REFERENCES 

Beskin V. S., 1990, Soviet Astron. Lett., 16, 286 

Beskin V. S., Gurevich A. V., Istomin Ya. N., 1993, Physics of the Pulsar Magnetosphere, Cambridge University Press 
Beskin V. S., 2005, " Osesimmetrichnye Stacionarnye Techeniya v Astrofizike" (Axisymmetric Stationary Flows in Astro- 
physics), Moscow, Fizmatlit (in Russian) 
Bisnovatyi-Kogan G. S., 1995, ApJS, 97, 185 

Bocquet, M., Bonazzola, S., Gourgoulhon, E., Novak, J., 1995, A & A, 301, 757 
Bonazzola S., Gourgoulhon E., 1996, A & A, 312, 675 
Boriakoff V., 1976, ApJ, 208, L43 

Cardall C, Y., Prakash M., Lattimer J. M., 2001, Apj, 554, 322 
Colaiuda A., Ferrari V., Gualtieri L., Pons J., 2008, MNRAS, 385, 2080 
Contopoulos I., Kazanas D., Fendt C, 1999, ApJ, 511, 351 

Courant R., Hilbert D., 1962, Methods of Mathematical Physics, John Wiley & Sons 
Ding K. Y., Cheng K. S, 1997, MNRAS, 287, 671 
Duncan R. C, 1998, ApJ, 498, L45 

Ginzburg V. L., Ozernoy L. M., 1964, Zh. Eksp. Teor. Fiz., 47, 1030 
Glampedakis K., Samuelsson L., Andersson N., 2006, MNRAS, 371, L74 
Goldreich P., Julian W. H., 1969. ApJ, 157, 869 
Gruzinov A., 2005, PRL, 94, 021101 

Hasskell B., Samuelsson L., Glampedakis K., Andersson N., 2008, MNRAS, 385, 531 
Israel G., et al, 2005, ApJ, 628, L53 
Komissarov S. S., 2006, MNRAS, 367, 19 

Landau L. D., Lifshitz E. M., 1987, The Classical Theory of Fields, 4th ed., Pergamon Press, Oxford 
Levin Y., 2007, MNRAS, 377, 159 
McKinney J. C, 2006, MNRAS, 368, L30 
Michel, F. C, 1973, ApJ, 180, L133 

Michel F. C, 1991, Theory of Neutron Star Magnetospheres, The University of Chicago Press 
Mofiz U. A., Ahmedov B. J., 2000, ApJ, 542, 484 

Morozova V. S., Ahmedov B. J., Kagramanova V. O, 2008, ApJ, in press (|arXiv:0806.2376|l 
Muslimov A. G., Tsygan A. I., 1992, MNRAS, 255, 61 
Pacini F., Ruderman M., 1974, Nat., 251, 399 
Pons J. A., Geppert U., 2007, A&A, 470, 303 
Rezzolla L., Ahmedov B. J., 2004, MNRAS, 352, 1161 

Rezzolla L., Ahmedov B. J., Miller J. C, 2001a, MNRAS, 322, 723; Erratum 338, 816 (2003) 



18 E. B. Abdikamalov, B. J. Ahmedov and J. C. Miller 



Ruderman M., Sutherland P. G., 1975, ApJ, 196, 51 

Samuelsson L., Andersson N., 2007, MNRAS, 374, 256 

Schwartz S. J., Zane S., Wilson, R. J., et al, 2005, ApJ, 627, L129 

Shapiro, S. L., Teukolsky, S. A., 1983, Black holes, White Dwarfs and Neutron Stars. Wiley, New York 

Sotani H., Kokkotas K. D., Stergioulas N., 2007, MNRAS, 375, 261 

Sotani H., Kokkotas K. D., Stergioulas N., 2008, MNRAS, 385, L5 

Spitkovsky A., 2006, ApJ, 648, L51 

Strohmayer T. E., 2008, AIPC, 968, 85S 

Strohmayer T. E., Watts A. L., 2006, ApJ, 653, 593 

Sturrock P., A., 1971, ApJ, 164, 529 

Thompson C, Duncan R. C, 1995, MNRAS, v275, 255 

Thompson C, Duncan R., C, 2001, ApJ, 561, 980 

Timokhin A. N., Bisnovatyi-Kogan G. S., Spruit H. C, 2000, MNRAS, 316, 734 
Timokhin A. N., 2006, MNRAS, 368, 1055 
Timokhin A. N., 2007, Ap & SS, 308, 345 

Timokhin A. N., Eichlers D., Lyubarsky Y., 2007, ApJ, 680, 1398 
Tsygan A. I., 1975, A&A, 44, 21 

Unno W., Osaki Y., Ando H., Sato H., Shibahashi H., 1989, Nonradial Oscillations of Stars, 2nd ed., University of Tokyo 
Press 

Watts A. L., Strohmayer T. E., 2006, ApJ, 637, L117 
Watts A. L., Strohmayer T. E., 2007, AdSR, 40, 10, 1446 

Woods P. M., Thompson C, 2006, in Lewin W. H. G., van der Klis M., eds, Compact Stellar X-ray Sources, Cambridge 
Univ. Press, Cambridge 



APPENDIX A: CALCULATION OF SS 



The derivative with respect to (f> of the right hand side of equation (|64|) must be equal to the derivative with respect to 9 of 
the left-hand side of (|65[) . Using this condition and expression (|49[) for S, we obtain 



Aa(dtSS)\ r =R = - 



Sti fin 

N6v r d r AnSo + d 9 A n S + — ^-r9 A n 5 o 

r r sin o 



+ (deisinOSv®) + d<t,5v^\ + Nd r deSo deSv r + N 9 d r dsSo d<t,Sv r 

rsmfl V / sm 6 



We now introduce a new function for shorthand: 



« W - -m Tt 



Using this notation, one can rewrite equations (|57|) and (|55|l in the following way 



6S em (r,t)= T ^fi-8s em (t), 



SS(t,r,6,4>) = Y, E 



M 3 



5sem(t)Y ern (8,<j>). 



1=0 m=-t 

Substituting the right hand side of the last equation into the left hand side of ()A1|I , we get 

2 °° 



JpY, E Qe(r) d t Ss em (t)£(e+l)Y em \ r=R = 



1=0 m—-l 



N5v f d r A n So + —deAnSo + -^H— A n S o 
r r sm o 



+ (d e (smeSv § ) + d+6v*) + NdrdeSo d e 5v f + -JL-d r d<f,So d^Sv* 

r sm 9 \ J sm 6 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 
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Multiplying now both sides of this equation by Y£, m , [8 , <j>) , integrating it over solid angle 4ir and using the orthogonality 
condition for spherical harmonics f Yp*, m ,Ye m dQ, = <5 f / f <5 m / m , we obtain 



d t 5si m {t) 



M J 



<'(ni)rV(r) h, 



8v e Sv' 1 ' 
NSv r d r AaSo + —deAnSo + ——d^AnSo 
r r sm 8 



(A6) 



+ — t-t; (de(sm65v S ) + d<t,5v^\ + NdrdgSo 8g5v r + ^ „ d T d<t,So d<t,5v 
rsmd V / sm 2 8 



N 



Y£ m ,(8,<j>)dn\r=R. 



In the case of toroidal oscillations of a NS in the mode (£',m'), the oscillation velocity components are given by (|70l> . 
Substituting the right-hand sides of (|70[) and (I72|) into (|A6|) . we get the following expression for dt5st m (t) for a NS with a 
dipole magnetic field 



dt8sim{t) = 



1 Bof R M 3 fj R 



1(1 + 1) Rq t (R) 
x J \dgYz m (sm 9 cos \ ~ e ^ cos sin + i^d^Ytm sin 8 sin xj — 



(A7) 



Y l' m> 

sin 8 



-dQ 



where fj B = r] R e '"r*. Substituting this into (|A4|) . we obtain the expression for 5S for a toroidal oscillation mode (£',m') of 
a NS a with dipole magnetic field 

00 t 



ck6S{r,t) = Y, E 



1=0 m=-e 



BoRf n fi R r 2 q e (r) 
'.{£+!) R 2 qt(R) 



(A8) 



x J ^deYe m (sin 8 cos x — ^ c °s 8 sin x) + ie^d^Yi-m sin 6* sin xj 



If the magnetic dipole moment is aligned with the oscillation mode axis (x = 0), one can easily show that 
B Rf R M 3 fj R - 



dt5sem(t) = 
dtSS 



£(£ + 1) 
im' 



R 2 qe{R) ' 



BoRfnVii 



r 2 q t , (r) 



where 5 pi, is the Kronecker tensor. 



(A9) 
(A10) 



APPENDIX B: SOLUTION OF THE * sc EQUATION FOR x = 

Substituting the right hand side of equation (|A10j) and expressions (|70[) into (|T9|) . we obtain a boundary condition for the 
function ^ sc for toroidal oscillations of a NS, for the case \ — 0: 

m' 2 Y t , m , {■»,<!>)- 



*scl'=« = - B oRf R Vn 



cosi?9 tf y^ m /(i?, 



d-d 



(Bl) 



£'{£' + l) sintf 

where $ is the integration variable. Substituting the expression for dtSS given by ()A10|) into equation l|78fl. we obtain an 
equation for the SC potential, ^ sc , for a toroidal oscillation mode (£' ,m') of a NS having a dipole magnetic field aligned with 
the oscillation mode axis 

rV(r) 



2r 2 qi (r) cos 8 d r (* sc ) + d r [r'qi (r)] sin 8 do (* sc ) - + BoRf^dr [r'qi (r)] R2 ^)^ (8, cb) = 



(B2) 



This is a first-order partial differential equation. According to a well-known theorem from the theory of such equations, 
(|B2[) is equivalent to the follow ing system of first-order ordinary differential equations (see, for example, Chapter II of Volume 
II of lCourant fc HilbertJ (jl962j) for a thorough discussion) 



dr 



dO 



2r 2 qi(r) cos# d r [r 2 qi(r)]sin8 

The characteristics of this system are 
fo = t , 
<Pi = </> , 

tpa = \/-r 2 q\ (r) sin 8 , 



dcj> dt 

r2 ii'( r ) v„ ,m a\ ~ _ 



2 + 1) 



BoRf R V R 



dr[r 2 qi (r)] q e ,(r) Y e , m ,(8(r), 



qi(r) R 2 q e >(R) cos 8(r) 



■ dr , 



(B3) 



(B4) 
(B5) 
(B6) 

(B7) 
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where 9(r) depends on the variable r due to ()B6I) : 

e(r) = arcsin( ^ ) . (B8) 
V V~ r 9i ( r ) / 



The integral of equation (|B2[) is an arbitrary function of ipo, ip±, ip2 and <ps: 

r((p , i, tp 2 , (fi-i) = . (B9) 
Using this, one can obtain an expression for the general solution of equation (|B2[) : 

2£'(£' + l) 

where "I>i is an arbitrary function that has to be determined from the boundary condition (|B1[) . Using (|B1|I . we can obtain a 
boundary condition for $i: 

* / ^^rvtw^w*' ^-*- (BU) 

where ^sc l r =fl i s given by (|B1|I . Now using equations (|B10|l and (|B11[) . we can write a final expression for the general solution 
of equation (|B2|) in the following form 

., - /; ^ *^y> * + - 1^)".*,.] , (bid 

where r' is the integration variable. The function 6{r) under this integral must be substituted by (|B8|) and, after performing 
the integration, the function 02 should be replaced by (|B7[) , The unknown function $2 [\/~ r 2 q±(r) sin 8, (f>, tj is determined 
using boundary condition ()B1|| : 



[\/-r 2 gi(r) sin<9, </>,t] \ r=R = -B Rf R fj R J 



costf 9tfl£/ m /(tf,<£) 



+ sin. M - ^ 

In order to obtain $2 for arbitrary r, one has to express the right hand side of this equation in terms of sin#, and then replace 
sin# by y / [r 2 qi(r)]/[R 2 qi(R)] x sin#. We then have our analytical expression for the function for toroidal oscillations of 
a NS. 



APPENDIX C: ENERGY LOSSES FOR M' £ MODES 



In this Section we solve equation (IB2|I and calculate 60 together with the total energy losses for non-axisymmetric toroidal 
oscillation modes for small 9o- Taking the limit of small 8 in the boundary condition (|B1|) for ^ sc , we obtain the following 
approximate expression for ^sc\r=R 



*sc|r=B = -B Rf R 7j R 



1 - 



£'(£' + !) 



t" m' 



,(2) 



ra'A™ , 



+ 4 2) , 

1 I'm' 



2(m' + 2) (to' + 2)£'{t' + 1) 



□ m +2 



Using this result, one can show that the general solution of equation (|B2I) for small 6 has the following form 



-B Rf R fj R gi(r) 
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(CI) 



(C2) 



where we have introduced a new function ge(r) = r 2 qi(r) for simplicity of notation. Using (|44p and (|96jl . we obtain an 
expression for A£ : 



AS = -B Rf R fj R N R [ cos d g Y t , m ,d6 
Jo 



(C3) 
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Substituting the expansion of Ye/ m i for small 9 given by (|106|l into (|C3[) , we get the following expression for A£ for small 



AS = -B Rf R fj R N R 



a™ ,e m ' + ( j&l , - , m , 4;' , ^ * m '+ 2 

t'rn' 1 \ I'm' n/ / . <~»\ f'm' 



2(m' + 2)' £ ' m ' 

Now substituting (|C2|I into (|46|) . we get an expression for p sc (R,8, (f>): 

Psc (R,e,4>) = -^fi D ^'{R,M)e m ' , 

where we have introduced a new quantity D^ m i (R, M) for simplicity of notation: 



(04) 



(C5) 



D t , m ,(R,M) = %ri4A$l,(l + m') 
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(C6) 



+ 2/(/ + l)Rg 1 {R)g l ,{R)dl r g 1 {r) 



Substituting (|C4|) and (|C5[) into (|103[) . we obtain an algebraic equation for 6>o for m/0 toroidal oscillation modes, which has 
the solution 



1QN R fj 2 R \A$l,D e , m ,(R,M)\ 



fl 

■I n 



(C7) 



We can see from this expression that Oq is small for modes with m' < 3, in agreement with the estimate of TBS. Remarkably, 
#o does not depend on Bq as in the case of m' = 0. Substituting ()C4[) . (|C5[) and (|105[) into (| 104f) gives 

2 



W = p ," (B Rf R v R ) 2 \A$l, D t , m ,(R,M)\ 8 2m ' +2 . 
Replacing (9o in this equation by the right hand side of equation (|C7|I . we obtain 
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